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Abstract 

To study the behavior of the Kazakov-Migdal at large N the quenched momen- 
tum prescription with constraints for treating the large N limit of gauge theories is 
used. It is noted that it leads to a quartic dependence of an action on unitary matrix 
instead of a quadratic dependence discussed in previous considerations. Therefore 
the model is not exactly solvable in the weak coupling limit. An approximation pro- 
cedure for investigation of the model is outlined. In this approximation an indication 
to a phase transition for d < 4, 8 with (3 cr = 8 is obtained. 
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1 Introduction 



Recently Kazakov and Migdal [Q] have made an important progress in an attacking the old 
problem of an evaluation of the large N limit for gauge models . They proposed a lattice 
gauge model which is solvable in the large N limit under an assumption of translation 
invariance of a master field and has a nontrivial critical behaviour. Different aspects of 



the Kazakov and Migdal (KM) model were considered - Jl0|. A relation of this model 



to QCD and, in particular, the property of asymptotic freedom still should be clarified. 



In this note a standard quenched momentum prescription [11] - [19| for treating the 



large N limit is used for the KM model . It has been established [T2|, [TJj that in the case 
of gauge theories the quenching of the momentum must be accompanied by a constraint 
on the eigenvalues of the covariant derivative. It will be shown that this constraint leads 
to a quartic dependence of an action on unitary matrix instead of a quadratic dependence 
discussed in previous considerations. In the strong coupling limit one can expect that the 
constraint can be omit and one lefts with an action which is quadratic on the unitary 
matrix V. In this case we can identify the translation invariant master field with a 
reduced field without quenching, i.e. without constraints. The relation of the translation 
invariant master field with the reduced field without quenching was noted by Makeenko 
|7J. However one cannot ignore the constraint in all regions of the coupling constant. This 
drops a hint that the approximation of the constant master field does not fully describes 
an asymptotic of the KM model at the large N limit. This constraint makes the KM 
model not exactly solvable for the weak coupling. We calculate the integral over unitary 
matrix in the semiclassical approximation, i.e. at small coupling. By the analogy with 
the Gross- Witten [[HJ] and the Brezin-Gross |5tJ models it is natural to expect a phase 
transition for a theory with an action being a polynomial on the one link unitary , in 
particular, for the quartic action. In the framework of some approximation procedure 
scheme this phase transition gives rise to a phase transition for the KM model. 

The paper is organized as follows. In Section 2 the general quenched momentum 
prescription with constraints is applied to the KM model and a reduced action with quartic 
dependence on the gauge field is presented. Then, in Section 3 we integrate out the gauge 
fields in the semiclassical approximation. In Section 4 we use an approximation scheme 
to evaluate the remaining integrals over the quenched momentum and the eigenvalues of 
the scalar field. The concluding remarks are collected in the last Section. 



2 Quenched Momentum Prescription for the KM model 

The KM lattice gauge model is defined by the partition function 



X,fl 



Here the field takes values in the adjoint representation of the gauge group SU(N) 
and the link variable U^(x) is an element of the group, \x = 1, ...d. 

The well known procedure for treating the large N limit of N x N matrix models con- 
sists in application of a quenched momentum prescription. This prescription for treating 
the limit of infinite N of theories with a global, or local gauge, U(N) symmetry, both on 
the lattice and in the continuum, has been obtained more than ten years ego by Eguchi 
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and Kawai [11 



Bhanot, Heller and Neuberger [12], Parisi [[HJ and Gross and Kitazawa 
Generally speaking, according to this prescription a reduced 



(see also [re- 
model can be described as containing just one side (one space-time point) and one scalar 
field $ on this side, or d links attached to one side and d fields for vector field. To 
get the reduced action one should replace a matrix field $(x) with D(x)^D Jf (x), where 
D(x) = e iPMXfl , and p^ is the diagonal matrix with matrix elements p l ^,i = 1, ...N. Then, 
one gets the vacuum energy in the large N limit by integration the free energy of obtained 
action over p. It has been shown that quenched theory produces the standard Feynman 
diagrams for invariant Green functions in all orders in perturbation theory. 

In the case of gauge theories the quenching of the momentum must be accompanied 
by a constraint on the eigenvalues of the covariant derivative. Without this constraint one 
gets naive reduced model without quenching. In the case of the Wilson gauge theory the 
naive reduced model describes correctly the theory at large N only in the strong coupling 
regime. 

Let us apply the quenched momentum prescription to the KM gauge theory ([!]). If 
we were to follow the quenched momentum prescription we would replace $(x + /i) with 
D fJi ^(x)Dj l , where = e iPtJ - a , and p M is the standard diagonal matrix whose elements are 
\p l | < n /a, a is the lattice spacing. The quenched model will then have an action 



and the vacuum energy 



F{p) = In 



9 fj,>o 



(2) 



(3) 



The vacuum energy per unit volume at the large N limit is then obtained by integrating 
F(p^) over all values of p M 

F= lim / dp,F(p,). (4) 



The integration over p^ is normalized to unity for all values of N 



d r w/2a N fld n , 1 

d -nL /2 m^](^. 



(5) 



If we were used the standard Haar measure, dU^, then by a change of variables U^ 
we would eliminate from the action (0). Thus no quenching would have occurred, and 
we would recover the KM model in the translationally invariant master field approxima- 
tion, 



In 



f d^dUy exp(a d iV[- tr V ($) + (3 J2 tr $ u ^ u t\) 



(6) 



The same procedure being applied to the standard Wilson theory yields to the Eguchi- 
Kawai reduced model without quenching, which is known to be necessary to obtain the 
correct results in weak coupling. The exit from this problem is known, one has to intro- 



duce a gauge invariant constraint on the U^s 



The suitable constraint restricts the 



eigenvalues of U^D^ to be equal to D^. Note that in continuum case the analog of these 
constraints yield the correct coupling constant renormalization [|14[]. We are going to use 
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these constraints for the reduced model (0). This means that we have to take the measure 
dfi(U) in (^) to be 

dfi{U) = l[dU tl C{U tM Dj, (7) 
where dll^ is the Haar measure on SU(N) and 

C(U„ D„) = n / dV, A(D,) 5(U, - V^VjDj) (8) 

A(D M ) = nsin 2 (^^a). 

For the moment we omit the Faddeev-Popov determinant which should be taken into 
account if one wants to get a correct result in weak coupling. If we integrate out the U^, 
setting Ufj, = V^D^V^D^ we obtain 

exp(F(p)) = / d1>Y[dV,exp{a d Ni-trV ($) + /3 £ tr ( QV^VftV^l V})]}. (9) 

To find the vacuum energy at the large N limit one should integrate F(p) over dp p. In 
fact we will use the formula 

F= j dfi{p)F{p), (10) 



where ^ 

Mp) = II T^r d ex p(-^ 2a2 )- ( n ) 

* \ a > 

In comparison with equation @, where we have the integral from the exponent con- 
taining the quadratic dependence from the unitary matrix U, in (||) we have to integrate 
the exponent with the quartic dependence from the unitary matrix. 

3 Weak Coupling for the Reduced KM Model with 
Quenching 

The action in ([]) contains the quartic interaction between the unitary matrix V and the 
hermitian matrix $, and the model does not look the exactly solvable. It is clear that the 
integral I($,D) 

/($,£>) = j dVexp[f3tr($VDV^VD j V^] (12) 

depends only on eigenvalues of the hermitian matrix $, 0j . 

We are going to calculate the integral over V^'s in ([5]) in the semiclassical approxima- 
tion. To this end let us note that the corresponding classical equation has a form 

DV^V^V^ - V^V^V^V,DV^+ 

DW^V^V^ - V^V.DV^V.DX = (13) 
Here (<fi)ij = S^fa. The solutions of this equation are 

Vq = TP, (14) 
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where T is any diagonal unitary matrix and P is any N x N permutation matrix which, 
when applied to an N- vector if), gives 

Pijipj = if) P (i). (15) 

Thus, for each permutation P, one has a solution of equation (0). So, for each P, we 
start with the decomposition of V around TP, 

V» = V ^, (16) 

where 

u = i + i\ a i a - ~\ a \ p u P + o(e) (17) 

To get the complete answer in the semiclassical approximation, we have to sum over all 
saddle points 

z{4>) = n / dv^*™ ~ nE / dut? dN w+> D >» p > u) (is) 

„ J fi p around P 

The effective action S(<f>, D u , P,U) up to quadratic terms on £ has the form 

su> i ,D„p,U) = Y,<l>m+ 

2 tr ( A • O0p(A • O0p + tr ( A • i)<p P D^ P {\ ■ 
+ tr ( A • £)L> M (A • 04>pDl<f> P - 2 tr ( A • 0<MW • 0^ - 2 tr ( A • O'(0p) 2 (19) 

where 0p is a diagonal matrix diagtpp = (f>pa),(f>p(2), ■■■4>p{n)\ £ ± ' 4 "' are the components of 
(A ■ £) along of generators of U(N) (from this place for simplicity we deal with U(N)) 

(a + *W = (44 - (A *W = + 4$); ( 2 °) 

(the components corresponding to generators being diagonal matrices do not make a 
contribution in (pl|) ). Taking into account the formula 



tr[A(\-QB(\.0-AB(\-$(\-0] = -^[A i B j +A j B i -A i B i -A^ 

i<j 



(21) 



we get 



S(<p t , D„ PM)=Y, <Ph) + 2 E(^« " <PpU)) 2v »iM +>ij ) 2 + (r*) 2 ] (22) 

i<j 

where 

V^^iD^-D^D^-D^). (23) 



In the basis (|2Cf), the measure for U(N) takes the form 



1 N -> N 



= n^'^-' iJ 'exp[--iV((e + «) 2 + (C> ij m Il^e X p[--iVe 2 + rEfi) 2 ]. (24) 

i<j D t=l D i=l 
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Performing the gaussian integration over £ ±,y , we get 

1 




PW 3 a <W/3(0 F(i) - P(J) ) 2 P W - + f J • (25) 

Having in mind that the distribution of the external momentum pi is invariant under 
permutation i — > P(i) one can for any given P renumerate the momentum and therefore 
one can claim that all permutations give the same contribution, so the final answer in 
semiclassical approximation up to some normalization factor has a form 

^^S ^-k^ (26) 

Including the contribution of the factor (E^) in we mix the orders in the perturbation 
theory and the 0(g 2 ) term (g 2 = f3~ l ) is not correct. To get the correct 0(g 2 ) answer the 
two-loop corrections should be computed. Note that the expression ( |26"D is well defined 
for some equal eigenvalues of $, that is in accordance with the well defined integral ( |T0| ) 
over compact manifold, and only singular point for finite N is the point f3 = oo, i.e. one 
cannot neglect the contribution of the second term in the denominator of fl26|). Of course, 
for infinite N the integral (pi]) can have a singular points (compare with phase transitions 
for the Gross-Witten pi] and the Brezin-Gross |f20| models). A hint to a phase transition 
for N = oo gives a zero of the denominator. 

Let us to compare the integral ( |10|) with the Itzykson-Zuber integral [^|, where one 
integrates the exponent of a quadratic form on V. If we calculate the Itzykson-Zuber inte- 
gral using the semiclassical approximation, then the corresponding classical solutions have 
the same form, but the contributions around given P depend explicitly on P. Summing 
over P provides a compensation of the singularities coming from coinciding eigenvalues 
of matrix <f), and in this case there is no reason to mix the orders in the perturbation 
theory, and, moreover, the remarkable fact is that summing over P the leading terms 
of the semiclassical approximation gives the exact expression found by different methods 



22, 23 



Substituting (EJ]) in (|J) we get the vacuum energy as the large N limit of the following 
expression 



2/PJi-PU 



(Vol) j dpM j^-^^mn ^\^Z + (27) 



M kj" 11 vr! -rji 1 6/3 



The free energy (|27| ) can be rewritten using the replica trick ||24|| . The replica trick 
allows us to estimate the free energy of the quenched system as the analytic continuation 
of the free energy in a class of annealed systems with the action being the copy of n's 
actions (127T). We have 
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T= lim- ( f dpf[d<p a e Nad ^< + W^ 2 & 2 (r)Y{Y{ 



A* i<j " V^i ^3> TO 1 6/3/ 

+ const. (28) 
Equation (^7]) points that there is a phase transition for some (3 = m^f3 cr . The reason 



is the following. For (3 = oo one has not the second term in the denominator of ( E7[) and 
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we left with the answer which was expected in the naive continuous limit and which is 
unstable, and for (3 = we get a stable model. So, it is natural to expect that the model 
exhibits a phase transition. 

To estimate the order of (3 cr one can roughly estimate (53) as 

T = In / ^ e ^H+W£^A 2 (^niI , 71 A ~ r, (29) 



where < (0j — 4>j) > means the average with semi-circular distributions of eigenvalues 
given by u(<j>) = - m 2 eff (j) 2 , i.e. 

iV 2 N 
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(0 - , ) 2 w(0)Wa(0')^0' = = —2~- ( 30 ) 



A J m eff m eff 



In our case m 2 ^ = a d (ml — d(3) = a d m^(l — df3 ), /3 is a dimensionless constant. 
< V^j > means the average of (^3|) with the measure (|TTD, so 



< >= 2(1 - 1/y/e). 
Therefore, the average of the first term in the denominator is equal to 

2(1-1/^)^ 
1-dPo 

So, for /3o = 1/d the contribution of the first term becomes infinite and we can neglect 
the second term and we get unstable model what is in accordance with expected phase 
transition for f3 Q = 1/d. 

The region f3 > 1/d is forbidden, since the action becomes unbounded from below. 
Let us examine the region f3o < 1/d including the negative (3q. In this case the action 
is bounded from below, but still there is a reason for a phase transition due to zero of 
denominator in (p9|) and 

fc - J^Ts- (31) 

This phase transition occurs only for d < 4,8. 

As for numerical factors we would like to note that they may be changed by a factor 
1/J which takes into account the two-loop corrections, i.e. in ( pB[) must be j instead of 
g. If J is such that 0, 8J > d then we can expect a phase transition. For negative J the 
above estimation gives a phase transition without any restrictions on d. 



4 Concluding Remarks 

In this paper an application of the quenched procedure to the KM model has been consid- 
ered. It has been noted that using the quenched procedure with constraints one obtains 
the one side model with the quartic interaction. This model has been investigated in the 
semiclassical approximation and the phase transition at [3 CT = 8 has been found. 
It will be shown in forthcoming paper how it is possible to modify the initial 



model to get the reduced model with quenching which admits an analytical treating at 
the large iV limit using the technique of the orthogonal polynomials. 
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